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The chiral Luttinger model for the edges of the fractional quantum Hall effect is obtained as 
the low energy limit of the Chern-Simons theory for the two dimensional system. In particular we 
recover the Kac-Moody algebra for the creation and annihilation operators of the edge density waves 
and the bosonization formula for the electronic operator at the edge. 
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The fractional quantum Hall effect (FQHE) is a man- 
ifestation of the formation of an incompressible state in 
a two dimensional electron system at discrete, density 
dependent values of the magnetic field. The major step 
towards a theoretical understanding of this phenomenon 
was taken by Laughlin fij] who proposed a variational 
many body wavefunction to describe a correlated, incom- 
pressible electron liquid at filling factors v = l/<f> for odd 
<j>. The microscopic wavefunction approach pLpj was sub- 
sequently augmented by an effective field theory of the 
Chern-Simons type. As it turned out both the bosonic 
|pf and fermionic B versions of the theory proved to be 
fruitful in explaining different aspects of the effect. While 
the incompressibility of the FQHE implies a gap for the 
excitations in the bulk of the system the edges support 
gapless modes that control the low energy physics of the 
problem ||. Wen pioneered the use of the chiral Lut- 
tinger model to describe these edge excitations and in- 
dicated, among others |p|-[10|, various non Fermi liquid 
effects which arise from this description. While plausible 
arguments, supported by recent numerical calculations 
Jll| , have been put forward in favor of the applicability 
of the Luttinger model to the edges of the FQHE §-|,|||, 
a microscopic justification is still lacking. It is the pur- 
pose of the following Letter to supply such a justification 
for the fractions of the form v — \/4>. In doing so we 
reveal an intimate relation between the bosonic Chern- 
Simons theory for the two dimensional system and the 
one dimensional Luttinger model of its edges. The lat- 
ter is obtained as the low energy limit of the former. 
In particular we recover the Kac-Moody algebra for the 
creation and annihilation operators of the edge density 
waves and the bosonization formula that enables one to 
express the electronic operator at the edge in terms of 
them. 

The central idea of the Chern-Simons theory is the 
statistical transmutation of the electronic field operators 
i/>(r) via an operator valued transformation 



_ ~*A(r 



(r) = e 



</>(r), A(r) = |dV/(r-r>(r'), (1) 
where p(r) = ip'(r)ip(r) = <fi'(r)(f>(r). We will be con- 



cerned with transformations that lead to bosonic </>(r) 
and <?r (r). It can be easily checked using the Baker- 
Hausdorff formula that this is the case when /(r — r') = 
/(r' — r) + <f)ir for odd <fi. The Chern-Simons transfor- 
mation also introduces additional terms in the equation 
of motion of <^(r) which couple to it as vector and scalar 
potentials |13| 
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where F(r — r') = V r /(r — r') and v(r) is the velocity 
operator. These "statistical gauge fields" constitute the 
main motivation for taking step (Jl|). For appropriate 
electronic densities they cancel on the average the exter- 
nal fields thus removing the degeneracy of the ground 
state and allowing for perturbative treatment of the in- 
teractions. In the Lagrangian formalism the equivalent 
bosonic system is described by the Chern-Simons action 
||. Using the polar decomposition of the bosonic field 
<j) = ^fpe ld and the velocity fields 

h e 

v l = —d l 6 + — {(n+Ai), (4) 

m mc 

with A taken such as to produce a constant magnetic 
field along z, one obtains the hydrodynamic form of the 
Lagrangian density which is more convenient for our pur- 
pose, 
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Henceforth, unless otherwise specified, length is mea- 
sured in units of the magnetic length I = ^Jhc/eB, time 
in units of inverse cyclotron frequency lu c = eB /mc, en- 
ergy is normalized by fiuj c and the density by the Lan- 
dau level degeneracy po — l/2nl 2 . On deriving (0) sev- 
eral surface terms were neglected as they can be shown 



to vanish for the configurations that will be of interest 
to us. Minimizing this action results in the mean-field 
equations which are the starting point of our analysis 
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We will consider a semi infinite system defined by an 
infinitely high wall situated at x < and obeying peri- 
odic boundary conditions over length L in the y direction. 
Our plan is to obtain the edge excitations as the random 
phase approximation (RPA) modes of the theory, i.e., 
the eigenmodes of the above equations linearized around 
a static solution. This step was carried out in jlj] and 
its essentials are reproduced below for the sake of com- 
pleteness of presentation. Once the gapless modes are 
at hand we will use them to expand the deviations of 
the various fields from their average values, namely the 
static solution. The expansion coefficients will become, 
once quantized, the density operators of the Luttinger 
model. Their commutation relations and Hamiltonian 
will be deduced from the quadratic Lagrangian obtained 
by expanding <M) to second order in the deviations. Fi- 
nally by inverting the transformation (nh we will recover 
the bosonization formula. 

The translational symmetry of the system along the y 
direction is in conflict with the canonical choice 



/(r-r') = 0tg- 1 



y-y 



x ■ 



X' 
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and the "symmetric gauge" induced by F(r — r') = 
lim e ^ + </>zx (r-r')/(|r-r'| 2 +e). Note that /(r - r') 
is a multi-valued function and we pick its branch cut 
along the positive x — x' axis. At this point we take ad- 
vantage of the fact that by adding to /(r — r') a function 
g(r — r') satisfying g(r — r') = g(r' — r) we do not affect 
the statistics of </)(r). Such a change on the operatorial 
level corresponds to a regular gauge transformation in 
the Lagrangian formalism. We choose 



5(r 



lim 



y-y 



*o+ ^{y-y') 2 + e 



y/{y-y') 2 +e 
(8) 



which is both single-valued and symmetric under the 
exchange of r and r'. As a result we find that (/ + 
g)(r — r') = ir<fr[l — e(y — y')] with a branch cut along 
the positive x — x' axis and V r (/ + g) = [O,27T0e(x — 
x')S(y — y')] which in turn implies the "Landau gauge" 
for the statistical potentials 



a(r) 



0, -4> d 2 r'e(x-x')8(y-y')p(r') 



(9) 
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with e(x) = 0(x) — 1/2 where O(x) is the step function. 
We will find it more convenient, however, to obtain the 
static solution and the edge modes in a slightly different 
gauge where we replace e(x — x') in (g, |l^) by 0(x — x'). 
Utilizing this gauge we look for a static solution in the 
form 



e = Ky- p(K)t , p = p(x) . 



(11) 



Due to the infinite confining potential we set the density 
to zero at the wall. For the external vector potential we 
use the gauge A = [0, Bx] which gives a velocity field of 



the form v x 
ansatz into (16 



and v y = v y (x). Inserting the above 



(6d) one finds a set of coupled equations 
for p(x), v y (x) and ao(x) that was solved numerically. 
For a fixed value of K the value of p was determined by 
requiring that p approaches its bulk density p = po/4> 
far from the wall. Under this condition p is the energy 
which is needed in order to add a particle to the edge. 
Representative examples of the density and current den- 
sity profiles for solutions (|ll]) in the case of short range 
interactions, along with p(K), are shown in Fig. |l|. One 
finds a one parameter family of static solutions depending 
on K and differing by the density of particles at the edge. 
Changing the value of K translates the condensate </> in 
the x direction. Using v y (0) = K and the fact that the 
velocity falls to zero in the bulk we can integrate (|6a|) to 
find that K/<p is the excess charge per unit length along 
the edge relative to a step-like constant density profile. 

Next we proceed to obtain the edge modes as solutions 
to the RPA equations, namely the linearized version of 
(pa)-(Sd) around one of the solutions (nTj) 
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where p denotes the part of the linearized " quantum pres- 
sure" term, i.e. the second term on the r.h.s. of (pa), 
containing x derivatives. The derivation of the edge exci- 
tation is facilitated by the observation that the derivative 
of the static solution with respect to K constitute a so- 



lution to (12a- 1 2d). This is the translation mode (along 
the x direction) of the system and it describes the addi- 
tion or removal of charge from the vicinity of the wall. 



Apriory the amount of charge carried by this solution is 
arbitrary since, due to the linearity of (12a-12d), it is de- 
termined only up to a multiplicative constant. We will 
find, however, at the end of the discussion, that the pe- 
riodic boundary conditions impose integer values for this 
charge in units of 1/0. We thus interpret the translation 
mode as the state of a smeared quasi-particle along the 
edge. For a representative example of these solutions see 
Fig. |. 

The gapless branch is obtained by modulating the 
translation mode in the y direction. Concentrating on 
the long-wavelength limit we find to lowest order in k 



dp = d K p(x)e' liky - UJt) 



Sv x = — 



<t>p{%) 

8v y = d K v y (x)e-^ ky -^ 

Sa = d K a {x)e- i{ - ky ' wt) 



ud K v y {x) - kd K a {x) - uAe'^-^ 



(13) 
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The real and imaginary part of (fL3|) constitute , to first 



order in k, independent solutions of (12a-12d), provided 



the dispersion relation u> = w(fe) is taken to be 



J 



w(A) 



dp 

dK 



lim — 



dx'd K p(x') U(x - x', k) - U{x - x', 0) 



(14) 



where U(x,k) — f dyU(x,y)e lky . For short range inter- 
action the last term in (|14j) vanishes while for Coulomb 
interaction it is proportional to — ln(fc) minus the elec- 
trostatic charging energy of the edge [Q . 




FIG. 1. Density in units of p (solid lines) and current den- 
sity in units of pu) c l (dashed lines) for the case of short range 
interaction U — {2-kK 2 /m)8{v — r'). The highest curve in 
each set corresponds to K = 2. Consecutive curves differ by 
AK = — 1. The right inset depicts the density and current 
density profiles of the translation mode corresponding to the 
K — 1 solution while the left one displays p(K). For more de- 
tails, especially on the differences between the cases of short 
range and Coulomb interactions, see Ref. |14[ . 

The solutions of the time independent linear eigenvalue 
problem defined by (12a- 12d), i.e. the same equations 
but with time derivatives replaced by iu>, constitute a 
complete set which one can use to expand the deviations 
of the various fields from their static average configu- 
rations. We will, however, restrict the sum over eigen- 
modes to the gapless edge excitations found above. In 



making this approximation we neglect other modes such 
as bulk magneto-plasmons and magneto-rotons |L5 16 
and their scattering states off the wall. We also do not 
consider the effects of quasi-particles and quasi-holes in 
the bulk. All of these excitations involve an energy gap 
and are thus expected to be irrelevant in the low en- 
ergy limit of the theory. At this stage we transform back 
to the original gauge (P,[lO|). As a result 9 changes to 
[K — (f)/2 J dxp(x)]y — \jj, — (f>/2 J dxpv v (x)]t and the ex- 
pansion of the deviation fields in the edge modes reads 



2vr 



2tt0 



*p = ^*p(o)0*p+^5>(*) fl - 



pe 



-iky 



fe^O 



27TI 

Tp~ 



5v x — —— 2_] Sp(k)[ujd K v y — kd K a — uj]e lky 
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Ar„ = ^Sp(0)d K v y + ^Y, Sp(k)d K v y e- lky (15) 
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The factor 2n<p/L in ( p^| ) was introduced, in connec- 
tion with the fact that J dxd K p(x) — l/4>, in order to 
make the configurations carry a unit charge (in dimen- 
sionful units) as we are interested to describe an elec- 
tron at the edge. The expansion coefficient Sp(0) will 
become, upon quantization, the number operator that 
measures the number of excess electrons on the edge 
relative to the static ground state configuration around 
which we expand p7|. It is time independent since we 



are interested in the case where the amount of parti- 
cles is fixed. The rest of the expansion coefficients 6p(k) 
are the new dynamical degrees of freedom and will play 
the role of the creation and annihilation operators of the 
edge density waves. Since the fields are real they obey 
5p(k) = Sp*(—k). The resulting effective low energy the- 
ory is described by the quadratic Lagrangian which is 
obtained by plugging ( |lq ) into (||) and keeping terms to 
lowest order in k. Here and in the remaining part of the 
Letter we restore the units of dimensions. 



L m 



.^£JW)*p(-*) 



k>0 



(16) 



*+ dfI :S P (0)S P (0) - 2 -^J2v(k)Sp(k)5p(-k) , 

fc>0 



L dK 



where v(k) = to{k)/k. The last two terms of ( fig ) are (up 
to a minus sign) the Hamiltonian of the reduced prob- 
lem. They coincide with the expression for the Hamil- 
tonian of the chiral Luttinger model ||,[l7J ■ The first of 
these terms corresponds to the change of energy, rela- 
tive to the ground state energy, due to excess charge on 
the edge. The second term describes the energy of the 
gapless excitations. The effects of intra-edge interactions 
are included in this Hamiltonian through the coefficient 
dp/dK and the dispersion relation v(k). The commuta- 
tion relations of the density operators are inferred from 



the symplectic part of L^ 2 ' , that is the first term in (16) 
We find that the conjugate momentum to 8p(k > 0) is 
— (2nh(f>i/kL)Sp(—k) thus giving, as anticipated, the al- 
gebra of the chiral Luttinger model 



kL 

[Sp(-p),5p(k)\ = — -j5 ktP 
Iirtp 



for k > 



(17) 



Note the factor 1/0 in (fijh which is responsible for the 
non Fermi liquid behavior of the theory j|-|§] . This alge- 
bra should be supplemented by a ladder operator U that 
raises the fermion charge at the edge in 1/0 steps [ p7| . 
Both U and the number operator Sp(0) commute with 
the rest of the operators. 

Within the low energy approximation and to low- 
est order in k we find that the bosonic operators, 
in the gauge defined by (P||lO|), are given by 0(r) = 
\J p{x) exp [i{9 + 59)(y)]. This expression can be used 
together with Eqs. (0) and (pj) in order to invert the 
Chern-Simons transformation^) to obtain an approx- 
imate form for the fermionic operators. The result con- 
sists, presumably, of the low energy edge components 
of the exact ip(r). Finally we obtain a one dimensional 
operator by projecting the two dimensional ip( r ) on the 
edge modes through an integration over the x direction 
with a weight given by the profile of their wavefunctions 



d K p(x)/y / p(x). The emerging expression is the bosoniza- 
tion formula for the electronic operator of the chiral Lut- 
tinger model §,0 



il>(y) = L 1/2 &cp[icp(y)] 
2m 



(18) 
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In the calculation leading to (ITq) y independent terms 
in A(r) were ignored. These terms do not arise in the 
case of an infinite system if one assumes that the den- 
sity vanishes at infinity. They do occur, however, for 
the case of a periodic system which we consider here, 
and they have to be omitted in order to make i/>(r) obey 
the correct anti-commutation relations. It is interesting 
that (p(y) is composed, in equal parts, of contributions 
from the Chcrn- Simon phase A and from the dynamical 
phase 69 of the bosonic field. We conclude by noting that 
the requirement that (|l8|) should obey periodic boundary 
conditions imposes a discrete spectrum for the number 
operator Sp(0), in units of 1/0. This observation is in 
accord with our earlier identification of the translation 
mode as the state of a quasi-particle at the edge. 
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